
*** Note that the p-values based on randomization inference vary slightly as would be expected***

use "$clean/clean_all.dta", clear
		   
loc experiments " canvass_any_main_unres  days_unres_main doors_unres_main_wins "

	   
		   
preserve

clear all
eststo clear
estimates drop _all

set obs 10
qui gen x = 1
qui gen y = 1

loc columns = 0

foreach choice in `experiments' `experiments' {

    loc ++columns
    qui eststo col`columns': reg x y

}

restore


/* Statistics */

loc colnum = 1
loc colnames ""

foreach choice in `experiments' {

/* Statistics */

 xi:  reg `choice' treatment , vce(r)
	estadd loc thisstat7 = `e(N)': col`colnum'
    sigstar treatment, prec(3)
    estadd loc thisstat2 = "`r(bstar)'": col`colnum'
    estadd loc thisstat3 = "`r(sestar)'": col`colnum'
	

	randcmd((treatment) reg `choice' treatment , robust), treatvars(treatment) reps(10000)
	mat define t=e(RCoef)
	local p_t: di %9.4f `=t[1,6]' 
	di `p_t'
	estadd loc thisstat4 = "[p=`p_t']": col`colnum'


	
	qui sum `choice'  if treatment==0
    estadd scalar thisstat6 = r(mean): col`colnum'

	estadd local thisstat8 = "No": col`colnum'
	
    loc ++colnum
    loc colnames "`colnames' `"`: var la `choice''"'"

	
 xi:  reg `choice' treatment $contr , vce(r)
	estadd loc thisstat7 = `e(N)': col`colnum'
    sigstar treatment, prec(3)
    estadd loc thisstat2 = "`r(bstar)'": col`colnum'
    estadd loc thisstat3 = "`r(sestar)'": col`colnum'
	

	randcmd((treatment) reg `choice' treatment $contr , robust), treatvars(treatment) reps(10000)
	mat define t=e(RCoef)
	local p_t: di %9.4f `=t[1,6]' 
	di `p_t'
	estadd loc thisstat4 = "[p=`p_t']": col`colnum'


	
	
	qui sum `choice'  if  treatment==0
    estadd scalar thisstat6 = r(mean): col`colnum'

	estadd local thisstat8 = "Yes": col`colnum'
	
	
    loc ++colnum
    loc colnames "`colnames' `"`: var la `choice''"'"


	}
	

	

loc rowlabels " "   "Treatment" " " " " " " "\midrule Control mean" "Observations" "\hline Controls"

loc rowstats ""

forval i = 1/8 {
    loc rowstats "`rowstats' thisstat`i'"
}


esttab * using "$tables/Table_1.tex", replace cells(none) booktabs nonotes nonum compress alignment(c) ///
 nogap noobs nobaselevels label stats(`rowstats', labels("`rowlabels'") ///
 fmt( 3 3 3 3 3 3 3 0 0 0 0 0 0 3 3 3 3 3 3 3 0  3 3 3 3 3 3 3 3 3 3 3 3 0 ))  mlabels(none) ///
   mgroups( "Any canvassing" "Days canvassed" "Door knocked",  pattern(1 0 1 0 1 0) prefix(\multicolumn{@span}{c}{) suffix(}) span erepeat(\cmidrule(lr){@span}))
